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Dispersion  of  Leachate  in  Aquifers 


PURPOSE:  Contaminated  dredged  material  is  often  placed  in  confined  disposal  facilities 
(CDFs)  designed  and  operated  to  control  environmental  impacts  of  the  disposed  sediment.  A 
CDF  is  a  diked  enclosure  having  structures  that  retain  dredged  material  solids.  When  contami¬ 
nated  dredged  material  is  placed  in  a  CDF,  contaminants  may  be  mobilized  to  form  leachate  that 
may  be  transported  to  the  site  boundaries  by  seepage.  The  purpose  of  the  research  presented  here 
is  to  examine  the  components  of  steady-state  leachate  attenuation  in  aquifers  and  to  develop  pre¬ 
dictive  equations  of  the  attenuation  for  use  in  a  screening  tool  being  developed  for  the  upland 
testing  manual  (USACE  2003).  The  main  factors  affecting  leachate  transport  and  dilution 
through  the  saturated  zone  of  an  aquifer  are  evaluated  to  develop  a  guidance  procedure  to  assist 
in  decision  making  regarding  the  use  of  leachate  controls  in  the  CDF.  The  U.S.  Environmental 
Protection  Agency’s  (EPA’s)  MULTIMED  model  is  used  to  develop  predictive  equations  for  the 
effects  of  recharge  and  lateral  and  vertical  dispersion  processes  on  center-line  concentrations  in 
the  aquifer.  Results  show  that  the  effects  of  these  processes  can  be  predicted  independently. 
Relationships  were  developed  to  estimate  the  attenuation  factor  for  each  process.  An  equation  for 
center-line  leachate  concentration  using  attenuation  factors  was  developed  to  predict  peak 
leachate  exposure  for  decision  making. 

BACKGROUND:  The  main  reason  for  the  study  described  herein  is  to  develop  a  set  of  simple 
equations  that  can  be  used  in  a  procedure  to  screen  the  leachate  pathway  for  unacceptable  expo¬ 
sure  to  contaminants  as  presented  in  the  upland  testing  manual  (USACE  2003).  Section  404  of 
the  Clean  Water  Act  of  1972,  as  amended,  the  National  Environmental  Policy  Act  of  1969,  the 
U.S.  Army  Corps  of  Engineers  (USACE)  management  strategy  for  dredged  material  disposal 
(Francingues  et  al.  1985),  and  the  USACE/USEPA  technical  framework  for  evaluating  the  envi¬ 
ronmental  effects  of  dredged  material  management  alternatives  (USACE/USEPA  1992)  require 
the  evaluation  of  the  confined  disposal  alternative  for  dredged  material  to  include  groundwater 
impacts.  The  leachate  pathway  is  depicted  in  Figure  1 . 

Leachate  seeping  into  the  groundwater  from  dredged  material  placed  in  a  CDF  is  produced  by 
several  potential  sources:  gravity  drainage  of  the  original  pore  water  and  ponded  water,  vadose 
zone  moisture,  inflow  of  groundwater,  and  infiltration  of  rainwater  and  snowmelt.  Thus,  leachate 
generation  and  transport  in  and  out  of  the  CDF  depend  on  many  site-specific  and  sediment-spe¬ 
cific  disposal  factors.  Leachate  concentrations  percolating  from  the  CDF  are  attenuated  as  they 
are  transported  through  the  vadose  zone  (Schroeder  and  Aziz  2002).  When  contaminants  reach 
the  saturated  zone,  they  are  transported  due  to  the  hydraulic  gradient  in  this  zone.  In  the  saturated 
zone,  contaminant  concentrations  are  a  function  of  the  dispersion  characteristics  of  the  aquifer, 
the  relative  rates  of  leachate  percolation  and  groundwater  flow,  and  the  aquifer  recharge  rates. 
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Figure  1 .  Model  of  dredged  material  leachate  pathway 


The  concentration  of  leachate  to  which  a  receptor  is  exposed  is  further  impacted  by  diffusion  or 
mixing  as  the  leachate  is  transported  from  the  CDF  locale  to  the  receptor  through  the  coarse¬ 
grained  layers  of  an  aquifer.  In  effect,  the  contaminant  concentration  of  the  leachate  is  diluted  by 
the  groundwater  flow.  Attenuation  by  adsorption  to  organic  matter  and  interactions  with  fine¬ 
grained  materials  will  also  occur  in  the  aquifer,  but  the  effect  is  generally  small  due  to  low  con¬ 
centrations  of  organic  and  clayey  materials  in  the  main  regions  of  saturated  groundwater  flow. 
However,  attenuation  of  leachate  concentration  in  the  saturated  zone  is  seriously  impacted  by 
groundwater  recharge  through  the  vadose  zone.  Hence,  leachate  at  the  receptor  is  affected  by  the 
groundwater  flow  rate  and  recharge  rate,  the  dominant  hydraulic  processes  between  the  CDF  and 
the  receptor.  Areas  with  high  groundwater  velocities  provide  greater  dilution  of  the  leachate 
plume,  but  spread  the  leachate  plume  more  quickly.  These  and  other  important  parameters  are 
discussed  below. 

The  quantity  and  quality  of  the  leachate  percolating  out  of  a  CDF  through  the  vadose  zone  and 
eventually  into  the  saturated  zone  is  a  function  of  CDF  design  and  operation,  the  contaminant, 
and  the  properties  of  the  vadose  zone.  Leachate  flow  through  the  vadose  zone  can  be  modeled 
using  the  HELPQ  model  (Schroeder  and  Aziz  1999).  Leachate  quantity  increases  with  the  area  of 
the  CDF  and  decreases  with  increased  application  of  best  management  practices.  Dewatering  and 
consolidation  of  the  dredged  material  decrease  the  pressure  head  that  drives  drainage  through  the 
CDF  and  decrease  the  hydraulic  conductivity  of  the  dredged  material,  both  serving  to  decrease 
leachate  production. 

Liners  and  drains  are  the  primary  control  features  for  leachate.  Liners  can  greatly  restrict 
leachate  flow  rates  from  CDFs  and  act  to  divert  leachate  to  offsite  drains  that  collect  the  leachate 
and  route  it  to  a  treatment  facility.  These  control  measures  prevent  nearly  all  of  the  leachate  from 
reaching  any  of  the  receptors. 

INTRODUCTION:  Leachate  that  reaches  the  saturated  zone  is  subjected  to  immediate  mixing 
and  dilution  with  the  groundwater.  The  rate  of  dilution  immediately  below  the  CDF  and  down- 
gradient  depends  on  the  relative  rates  of  leachate  infiltration  into  the  groundwater  flow  and  aqui¬ 
fer  properties  such  as  hydraulic  conductivity,  thickness,  porosity,  and  dispersivity.  Additionally, 
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aquifer  recharge  due  to  rainfall  infiltration  provides  additional  leachate  dilution  downstream  of 
the  facility.  The  attenuation  of  leachate  in  the  aquifer  is  a  three-dimensional  solute  transport 
process  that  is  represented  by 


—  W 

*  dx2  y  By2 


where  x,  y,  and  z  are  the  longitudinal,  lateral,  and  vertical  directions,  respectively  [L];  C  is  the 
dissolved  concentration  of  the  contaminant  [ MIL3 ];  D*,  D,  ,  and  Dz  are  the  dispersion  coefficients 
in  the  x,  y  and  z  directions,  respectively  [L2/T\;  V  is  the  uniform  groundwater  velocity  in  the  x- 
direction  [1/7];  t  is  the  elapsed  time  [7];  Vr  is  the  net  recharge  percolating  directly  into  the  con¬ 
taminant  plume  [L/T];  B  is  the  thickness  of  the  saturated  zone  [7];  and  8  is  the  effective  porosity 
[dimensionless]. 

The  flow  domain  is  regarded  as  semi-infinite  in  the  x-direction  (0  <  x  <  <x>),  infinite  in  the  y- 
direction  (-00  <y  <  oo),  and  finite  in  the  z-direction  (0  <  z  <  B).  A  schematic  of  the  flow  domain 
in  Figure  2  shows  the  coordinate  system,  and  the  facility,  receptor,  and  variable  definitions. 

There  are  a  number  of  numerical  models  that  can  be  used  to  simulate  this  process,  such  as 
MEPAS  (Whelan  et  al.  1996),  MULTIMED  (Salhotra  et  al.  1993),  AT123D  (Yeh  1981),  and  the 
Department  of  Defense  Groundwater  Modeling  System  (GMS),  to  name  a  few.  MULTIMED  is 
used  as  the  computer  model  for  the  evaluation  of  the  impact  of  aquifer  properties  on  leachate 
concentrations  at  the  source  because  it  has  been  widely  used  for  similar  applications  and  com¬ 
pletely  describes  the  dominant  processes.  In  solving  Equation  1,  MULTIMED  uses  two  bound¬ 
ary  conditions  in  each  of  the  x,  y,  and  z  directions  and  an  initial  condition. 

The  source  is  defined  as  the  downgradient  edge  of  the  waste  disposal  unit  and  is  used  as  the 
upstream  boundary  condition.  The  MULTIMED  model  allows  a  choice  between  two  boundary 
conditions  with  respect  to  the  distribution  of  contaminant  along  the  vertical  plane  at  the  source 
(Salhotra  et  al.  1993).  The  first  boundary  condition  specifies  the  contaminant  concentration  as  a 
Gaussian  distribution  in  the  lateral  direction  and  uniform  over  the  source  penetration  depth  H. 

Mathematically,  this  boundary  condition  can  be  expressed  as 

' 

C0  exp 

C(0,y,z,t)  =  < 

0 

in  which  C0  [MIL3]  is  the  maximum  dissolved  concentration  of  the  solute  at  the  source  and 
occurs  at  the  center  of  the  Gaussian  distribution.  The  standard  deviation  a  is  a  measure  of  the 
source  width  W  normal  to  the  mean  flow  direction  and  is  written  as 
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a 


W_ 

6 


(3) 


which  implies  that  99.86  percent  of  area  under  the  Gaussian  source  falls  within  the  width  of  the 
facility. 


The  concentration  at  the  center  of  the  Gaussian  distribution  is  determined  based  on  the  dilution 
of  the  contaminant  concentration  entering  the  saturated  zone.  MUMLTIMED  uses  the  mass  bal¬ 
ance  of  contaminants  to  compute  the  source  concentration  C0  as  (Salhotra  et  al.  1993): 


Afvf 

(AfVf+  Jbr  V  OH  a) 


where  Af  \s  the  area  of  the  facility  [/.'] ;  Vj  is  the  percolation  rate  from  the  facility  \LIT\\  Cf  is  the 
concentration  of  leachate  percolating  from  vadose  zone  [M/L^]',  and  H  is  the  source  thickness  [T].- 

The  thickness  of  the  source  represents  the  depth  penetration  below  the  downstream  end  of  the 
facility  that  is  impacted  by  the  leachate.  The  depth  of  penetration  is  impacted  by  the  vertical 
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advection  of  water  as  it  moves  beneath  the  facility  and  by  the  vertical  dispersion  in  the  saturated 
zone.  This  thickness  is  given  by  Salhotra  et  al.  (1993)  as 


H  =  (2  azL)112  +  B 


1  -  exp 


iV 

V9B 


(5) 


where  a2  is  the  vertical  dispersivity  [L]  and  L  is  the  length  of  the  facility  [I], 

The  first  term  on  the  right-hand  side  of  Equation  5  represents  the  contribution  of  vertical  disper¬ 
sion  to  the  depth  of  penetration,  and  the  remaining  terms  represent  the  effect  of  advection. 
Clearly,  if  the  source  thickness  H  exceeds  the  thickness  of  the  saturated  zone  B,  then  H  is  set 
equal  to  B.  In  this  case,  the  entire  saturated  zone  downgradient  of  the  facility  is  completely 
mixed  and  vertical  dispersion  becomes  insignificant  in  the  transport  process  along  the  vertical 
plane  of  symmetry. 

The  steady-state  solution  for  the  semi-infinite  domain  assumes  zero  concentration  at  an  infinite 
distance  from  the  source,  i.e.  the  background  concentrations  are  zero. 

C(oo,y,z)  =  0  (6) 

C(x,  ±oo,  z)  =  0  (7) 

The  final  boundary  conditions  assume  that  there  is  no  contaminant  flux  across  the  upper  and 
lower  boundaries  of  the  saturated  zone.  These  boundary  conditions  are  represented  as 

^r.(x>y>  °)  =  0  (g) 

dz 


^-{x,y,B)  =  0 

dz 

The  steady-state  concentration  along  the  plume  centerline  can  be  expressed  as 


(9) 


J_J  CO 

C(x,0,0)  =  £  —  J  exp 
B  0 


2  2 
a  u 


(u2  D 


-  x 


D, 


\l/2 


4D, 


du 


(10) 


where 


£  = 


exp 


f  Vx  ^ 
v2D*y 


M/2 

and  u  is  the  integration  variable  [Z/1]. 


(11) 
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Assumptions  used  in  the  analysis  are  as  follows: 

1.  A  single  aquifer  with  uniform  thickness  is  modeled.  The  saturated  porous  medium 
properties  are  isotropic  and  homogeneous. 

2.  The  groundwater  flow  velocity  is  steady  and  uniform.  This  implies  that  the  recharge 
through  the  facility  and  into  the  groundwater  plume  is  small  compared  to  the  natural 
(regional)  flow. 

3.  No  contaminant  degradation/transformation  or  sorption  is  assumed. 

4.  The  background  contaminant  concentration  in  the  aquifer  is  zero. 

AQUIFER  IMPACT  ON  LEACHATE.  The  screening  procedure  is  developed  based  on  the 
CDF  size,  saturated  flow  properties,  aquifer  properties,  recharge  rates  and  receptor  location.  The 
main  concern  of  this  screening  procedure  is  to  determine  the  peak  contaminant  concentration 
reaching  the  receptor  under  steady-state  conditions.  Due  to  the  symmetrical  nature  of  the  prob¬ 
lem,  peak  concentrations  will  occur  in  the  longitudinal  direction  along  the  centerline  of  the  facil¬ 
ity.  Therefore,  a  receptor  located  downstream  of  the  facility  along  the  centerline  of  the  facility 
will  experience  a  higher  concentration  than  all  points  that  are  at  the  same  distance  away  but  are 
not  along  the  centerline.  This  concentration  is  represented  as  a  function  of  several  variables. 


dh 


C  =  f"  C0  D,L,W,B,K,—,Vf,Vr,ay,az 

dxJ 


(12) 


where  D  is  the  distance  from  facility  to  receptor  along  centerline  of  facility  [L];  K  is  the  hydrau¬ 
lic  conductivity  [1/7];  dhldx  is  the  slope  of  the  water  table  [L/L];  and  ay  is  the  transverse  disper- 
sivity  [L], 

Equation  12  may  be  written  in  dimensionless  form  as 


f 


W 


B 


Qr 


Jcc^D  ’  ^azL  ’  L+D  qs+qf  qs+(lr 


(13) 


where  qs  equals  BK{dh/dx),  which  is  the  groundwater  flow  per  unit  width  of  aquifer  [1/7];  qf 
equals  V/L,  which  is  the  leachate  flow  rate  downward  per  unit  width  of  facility  [L  /7];  and  qr 
equals  VrD,  which  is  the  groundwater  recharge  rate  per  unit  width  [L2/T\. 

Since  C0  as  given  by  Equation  4  can  be  written  as  a  function  of  qf/{qs  +  qf),  and  H  as  given  in 
Equation  5  is  a  function  of  az,  Equation  13  can  be  reduced  further  to 


C  f  W  H_  L  qs 
C0  ^Ja~D  ’  B  ’  L  +  D  qs  +  qr  ^ 


(14) 


The  source  thickness  H  refers  to  the  maximum  depth  of  mixing  at  the  downstream  edge  of  the 
facility  and  is  defined  by  Equation  5.  There  are  three  plausible  scenarios  for  the  depth  of  pene¬ 
tration  in  relation  to  the  facility  and  the  receptor  as  depicted  in  Figure  3. 
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In  the  scenario  labeled  1  in  Figure  3,  the  source  thickness  is  greater  than  the  thickness  of  the 
saturated  zone  implying  that  complete  vertical  mixing  occurred  beneath  the  facility.  In  this  case, 
the  solution  assumes  H-B ,  and  the  plume  covers  the  entire  depth  of  the  saturated  layer  from  the 
downgradient  boundary  of  the  facility  to  the  receptor.  The  contribution  of  groundwater  recharge 
will  simply  have  a  dilution  impact  on  the  plume. 

The  second  scenario,  represented  by  the  vertical  penetration  depth  of  the  plume  (line  2  in  Fig¬ 
ure  3),  indicates  that  complete  vertical  mixing  of  the  saturated  zone  does  not  occur  beneath  the 
facility.  However,  complete  vertical  mixing  occurs  upstream  of  the  receptor.  Vertical  dispersion 
and  groundwater  recharge  serve  to  force  the  plume  to  become  completely  mixed  throughout  the 
depth  of  the  saturated  zone  before  it  reaches  the  receptor.  The  solution  of  this  scenario  is  the 
same  as  that  of  case  1 ,  since  the  entire  depth  of  groundwater  flow  is  enveloped  in  the  plume. 

The  third  and  last  scenario  is  the  case  where  complete  vertical  mixing  of  the  saturated  zone 
occurs  downstream  of  the  receptor.  This  is  a  case  of  one  or  more  of  the  following:  thick  aquifer, 
short  distance  to  the  receptor,  low  recharge  rate,  or  small  vertical  dispersivity.  In  this  case,  all 
parameters  of  Equation  14  are  significant.  The  thickness  of  the  plume  at  the  receptor  H’  is 
defined  as 


H'={2az{L+D))V2  +B 


f 

1  -  exp 

v 


LVf+DV^ 
V0B  , 


(15) 


The  actual  relationship  described  by  Equation  14  is  derived  using  MULTIMED  simulations.  The 
impact  of  aquifer  properties  on  the  leachate  concentration  at  the  receptor  is  evaluated  using 
MULTIMED  to  develop  a  screening  procedure.  The  model  could  be  run  for  the  site-specific 
conditions  if  increased  accuracy  in  the  predictions  were  needed  to  pass  the  screening. 

In  order  to  determine  steady-state  contaminant  concentrations  at  the  receptor,  MULTIMED  was 
used  to  simulate  the  effects  of  the  parameters  of  Equation  14  on  the  leachate  concentration  at  the 
receptor.  Simulation  parameters  may  be  divided  into  facility  (L,  W,  Q  qj),  aquifer  ( B ,  K,  dh/dx, 

ay,  az,  qs),  and  recharge  (Vr  and  D).  Table  1  summarizes  data  used  in  the  simulation.  The 
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selection  of  values  for  the  simulations  was  designed  to  isolate  the  effects  of  the  attenuation  proc¬ 
esses  as  explained  in  subsequent  sections. 


Table  1 

MULTIMED  Simulation  Data 


a2  B 
(m)  (m) 


1 

0.0252 

0.2 

500 

10 

31500 

0.003 

100 

1.89 

1.55 

15 

2 

0.0252 

0.2 

500 

20 

31500 

0.003 

100 

1.89 

1.55 

15 

3 

0.0252 

0.2 

500 

50 

31500 

0.003 

100 

1.89 

1.55 

15 

4 

0.0252 

0.2 

500 

75 

31500 

mm 

100 

15 

5 

0.0252 

0.2 

31500 

om 

100 

1.89 

1.55 

15 

6 

0.0252 

0.2 

500 

250 

31500 

0.003 

100 

1.89 

1.55 

15 

7 

0.0252 

0.2 

31500 

0.003 

100 

1.89 

1.55 

15 

8 

0.0252 

0.2  i 

31500 

0.003 

100 

1.89 

1.55 

15 

0.003 


0.003 


0.003 


0.003 


0.003 


3 _ 

3 


0.003 


0.003 


0.003 


0.0001 


160000 

0.00025 

106667 

0.0005 

66667 

0.001 

31511 

0.003 

21333 

0.005 

12000 

0.01 

2667 

0.05 

0.1 

i 

1.55  15 


1.55  80 


1.55  70 


1.55  50 


1.55  40 


1.55  35 


1.55  30 


1.55  60 


1.55  90 


1.55  50 


1.55  15 


1.55  15 


1.55  15 


1.55  15 


1.55  15 


1.55  15 


1.55  15 


1.55  15 


1.55  15 


0.01 

1.55 

6 

0.01 

1.55 

11 

0.01 

1.55 

14 

0.01 

1.55 

15 

0.01 

1.55 

20 

0.01 

1.55 

16 

0.01 

1.55 

28 

0.01 

1.55 

35 

0.01 

1.55 

42 
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Lateral  Dispersion.  The  rate  at  which  the  contaminant  spreads  in  the  lateral  direction  is  a 
function  of  the  lateral  dispersivity.  The  magnitude  of  lateral  dispersion  is  a  function  of  the  dis¬ 
tance  that  the  plume  travels,  and  the  reduction  of  centerline  concentration  due  to  lateral  disper¬ 
sion  is  a  function  of  the  facility  width.  The  results  of  MULT1MED  simulation  runs  1  through  9 
with  constant  recharge  and  vertical  dispersion  are  shown  in  Figure  4  to  evaluate  the  impact  of 
lateral  dispersion.  Figure  4  clearly  indicates  that  as  the  width  of  the  facility  increases,  with  all 
other  variables  remaining  the  same,  the  concentration  at  the  receptor  increases  and  approaches 
the  source  concentration  as  W  approaches  infinity.  Similarly,  for  narrow  facilities  (small  W)  the 
concentration  at  the  receptor  is  smaller.  Figure  4  also  indicates  that  the  increase  in  D  has  the 
opposite  effect  as  the  increase  in  W.  Lateral  dispersion  has  little  effect  on  peak  concentrations  of 
leachate  at  the  receptor  location  except  where  the  distance  to  the  receptor  is  very  long  (D  >  10 
W)  or  where  the  width  of  the  disposal  facility  is  very  narrow  (W  <  10  m).  The  lateral  dispersion 
factor  Fy  shown  as  a  solid  line  in  Figure  4,  was  determined  by  regression  to  be 


0.99562  +  8.95506 

l 


w 


-1.5 


Figure  4.  Impact  of  facility  width  on  centerline  concentration  (simulation  runs  1-9) 


Vertical  Dispersion.  As  explained  previously,  there  are  three  possible  cases  that  should  be 
considered  in  assessing  the  impact  of  vertical  dispersion  on  contaminant  concentration  at  the 
receptor.  These  cases  are  shown  in  Figure  3  and  illustrate  the  relationship  among  vertical  disper¬ 
sivity,  source  thickness,  aquifer  thickness,  and  receptor  distance  from  facility.  The  first  case 
represents  the  condition  where  the  entire  depth  of  the  saturated  zone  is  completely  mixed  verti¬ 
cally  beneath  the  facility.  In  this  case  it  is  assumed  that  H  =  B.  Since  C0  is  inversely  proportional 
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to  B,  the  concentration  at  the  receptor  is  inversely  proportional  to  B.  However,  the  dilution  by 
recharge  decreases  with  increasing  aquifer  thickness.  In  this  case  the  vertical  dispersion  is  so 
large  under  the  facility  that  it  has  no  impact  on  the  leachate  concentration  downgradient  of  the 
facility.  Hence,  the  vertical  dispersion  attenuation  factor/^  is  equal  to  1. 

The  second  case  shown  in  Figure  3  is  for  H  <  B  <  H’  where  H’  is  the  unrestricted  thickness  of 
the  plume  at  the  receptor  as  calculated  by  Equation  15.  In  this  case,  complete  vertical  mixing  in 
the  entire  saturated  zone  takes  place  at  a  point  between  the  downstream  boundary  of  the  facility 
and  the  receptor.  Vertical  dispersion  as  well  as  dilution  of  the  plume  by  the  recharge  contribute 
to  the  concentration  at  the  receptor.  The  concentration  at  the  receptor  is  inversely  proportional  to 
the  aquifer  thickness  and  is  independent  of  the  vertical  dispersivity.  To  account  for  complete 
mixing  occurring  between  the  facility  and  the  receptor,  a  vertical  dispersion  attenuation  factor  Fz 
can  be  applied  to  compute  the  concentration  at  the  receptor. 

Fz=-  07) 

z  B 


where  H  is  calculated  using  Equation  5. 

The  third  case  shown  in  Figure  3  indicates  that  complete  vertical  mixing  is  not  accomplished 
between  the  facility  and  the  receptor.  In  other  words,  H  <  H’  <  B.  Hence,  both  vertical  disper¬ 
sion  and  recharge  are  significant  factors  in  determining  the  concentration  at  the  receptor.  The 
concentration  at  the  receptor  is  not  impacted  by  the  aquifer  thickness,  but  is  a  function  of  the 
length  of  the  facility  and  the  distance  between  the  receptor  and  the  facility.  In  this  case  the  verti¬ 
cal  dispersion  is  not  restricted  by  the  aquifer  thickness  and  the  impact  of  vertical  dispersion  as  a 
function  of  receptor  distance  and  source  length.  The  results  of  MULTIMED  simulation  runs  10 
through  18  without  recharge  and  with  constant  lateral  dispersion  are  shown  in  Figure  5  to  evalu¬ 
ate  the  impact  of  vertical  dispersion.  The  vertical  dispersion  attenuation  factor  Fz,  shown  as  a 
solid  line  in  Figure  5,  was  determined  by  regression  to  be 
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Figure  5.  Impact  of  vertical  dispersivity  for  very  thick  aquifers  on  centerline  concentrations 
(simulation  runs  10-18) 


Impact  of  Recharge.  Groundwater  recharge  from  rainfall  infiltration  serves  to  dilute  the 
plume  and  hence  reduces  the  concentration  of  the  contaminants  reaching  the  receptor.  The  results 
of  MULTIMED  simulation  runs  19  through  27  with  constant  lateral  and  vertical  dispersion  are 
presented  in  Figure  6  to  show  the  effect  of  recharge  rates  on  center-line  concentrations.  The 
results  clearly  show  that  the  increase  in  recharge  rates  decreases  concentrations  at  the  receptor. 
Aquifer  recharge  has  little  effect  on  peak  concentrations  of  leachate  at  the  receptor  location 
except  where  the  distance  to  the  receptor  is  very  long  ( D  >  500  B ),  where  the  thickness  of  the 
aquifer  is  very  shallow  (B  <  1  m),  or  where  the  groundwater  water  velocity  is  very  slow  (V  < 
1  m/yr).  The  recharge  attenuation  factor  Fr,  shown  as  a  solid  line  in  Figure  6,  was  determined  by 
regression  to  be 


Fr  =1-0.9877 


Or 

qs(H'!B)+qr 


(19) 


Combined  Effects.  By  combining  the  results  of  the  analysis  presented  above,  a  relationship 
between  center-line  concentrations  at  the  receptor  can  be  represented  as 

p :-  =  FyFzFr  (20) 

and  is  represented  graphically  in  Figure  7  using  MULTIMED  simulation  runs  1  through  36.  Note 
that  runs  28  through  36  represent  simulations  using  other  values  of  lateral  and  vertical  dispersion 
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in  combination  to  provide  a  verification  of  the  overall  approach.  Equation  20  or  Figure  7  can  be 
used  to  predict  the  maximum  concentrations  expected  at  the  receptor. 


qr/[qs  (H’/B)  +  qr] 


Figure  6.  Impact  of  recharge  on  center-line  concentrations  (Simulation  runs  19-27) 


Figure  7.  Combined  effect  of  dispersion  and  recharge  (Simulation  runs  1-36) 
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Relating  C  to  the  concentration  of  leachate  entering  the  saturated  zone  C/,  Equation  4  is  substi¬ 
tuted  into  Equation  20  to  obtain 


C  = 


A/Vf 


( Af  Vf  +  n  VOHtJ 


'f 


Fy  Fz  Fr 


(21) 


where  B  is  substituted  for  H  when  H>  B.  Equation  21  is  an  equation  that  can  be  readily  solved  in 
a  screening  procedure.  It  requires  the  physical  parameters  of  the  aquifer  and  the  CDF  and  the 
leachate  concentration  entering  the  aquifer,  which  can  be  obtained  from  a  CDF  leachate  genera¬ 
tion  model  such  as  HELPQ  (Schroeder  and  Aziz  1999). 

CONCLUSIONS:  The  research  presented  above  provides  predictive  equations  of  leachate 
attenuation  in  aquifers  for  use  in  the  leachate  screening  procedure  presented  in  the  upland  testing 
manual  (USACE  2003).  Key  components  of  an  aquifer  such  as  thickness,  recharge  rate,  and 
aquifer  flow  rate  and  facility-related  parameters  such  as  infiltration  rate,  contaminant  concentra¬ 
tion,  and  facility  size  are  used  to  determine  the  maximum  contaminant  concentration  at  a  recep¬ 
tor.  The  procedure  accounts  for  dilution  of  the  leachate  under  the  disposal  facility  and  downgra- 
dient  dilution  by  recharge  and  lateral  and  vertical  dispersion.  The  relationships  developed  in  this 
technical  note  can  be  used  to  estimate  the  peak  concentrations  of  leachate  that  would  reach  the 
receptor  under  steady-state  conditions. 

Peak  concentrations  of  leachate  at  the  receptor  location  are  strongly  affected  by  vertical  disper¬ 
sion,  impacting  the  dilution  of  the  leachate  under  the  disposal  facility  and  downstream  of  the 
facility  until  the  leachate  plume  is  fully  mixed  throughout  the  entire  thickness  of  the  aquifer.  The 
attenuation  is  greatest  in  thick  aquifers  with  large  groundwater  velocities.  Lateral  dispersion  has 
little  effect  on  peak  concentrations  of  leachate  at  the  receptor  location  except  where  the  distance 
to  the  receptor  is  very  long  or  where  the  width  of  the  disposal  facility  is  very  narrow.  Similarly, 
peak  concentrations  of  leachate  are  largely  unaffected  by  aquifer  recharge  except  where  the  dis¬ 
tance  to  the  receptor  is  very  long,  where  the  thickness  of  the  aquifer  is  very  shallow,  or  where  the 
groundwater  water  velocity  is  very  slow. 

POINTS  OF  CONTACT:  For  additional  information,  contact  Dr.  Paul  R.  Schroeder,  (601-634- 
3709,  Paul.R.Schroeder@erdc.usace.army.miP)  or  the  Program  Manager  of  the  Dredging  Opera¬ 
tions  and  Environmental  Research  Program,  Dr.  Robert  M.  Engler,  (601-634-3624, 
Robert. M.Engler@erdc.usace. army.mil).  This  technical  note  should  be  cited  as  follows: 

Schroeder,  P.  R.  and  Aziz,  N.  M.  (2004).  “Dispersion  of  leachate  in  aquifers,”  DOER 
Technical  Notes  Collection  (ERDC  TN-DOER-C34),  U.S.  Army  Engineer  Research  and 
Development  Center,  Vicksburg,  MS.  www.wes.army.mil/el/dots/doer 

REFERENCES 

Schroeder,  P.  R.  and  Aziz,  N.  M.  (1999).  “Contaminant  leaching  model  for  dredged  material  disposal  facilities,” 
ASCE  Journal  of  Environmental  Engineering,  125  (9),  835-844. 

Schroeder,  P.  R.,  and  Aziz,  N.  M.  (2002).  ’’Effect  of  vadose  zone  on  CDF  leachate  concentrations,”  ASCE  Dredging 
2002,  Orlando,  FL. 


13 


ERDC  TN-DOER-C34 
July  2004 

Francingues,  N.  R.(  Jr.,  Palermo,  M.  R.,  Lee,  C.  R.,  and  Peddicord,  R.  K.  (1985).  ’’Management  strategy  for  disposal 
of  dredged  material:  contaminant  testing  and  controls,”  Miscellaneous  Paper  D-85-1,  U.S.  Army  Engineer 
Waterways  Experiment  Station,  Vicksburg,  MS. 

Salhotra,  A.  M.,  Mineart,  P,  Sharp-Hansen,  S.,  Allison,  T„  Johns,  R.,  and  Mills,  W.  (1993).  “Multimedia  exposure 
assessment  model  (MULT1MED)  for  evaluating  the  land  disposal  of  wastes,  Volume  1:  Model  theory,”  Rep. 
prepared  for  the  Environmental  Research  Laboratory.  U.S.  Environmental  Protection  Agency,  Athens,  GA. 

U.S.  Army  Corps  of  Engineers.  (2003).  “Evaluation  of  dredged  material  proposed  for  disposal  at  island,  nearshore, 
or  upland  confined  disposal  facilities  —  Testing  manual,’  Technical  Report  ERDC/EL  TR-03-1,  U.S.  Army 
Engineer  Research  and  Development  Center,  Vicksburg,  MS. 

U.S.  Army  Corps  of  Engineers/U.S.  Environmental  Protection  Agency  (USACE/USEPA).  (1992).  ’’Evaluating  envi¬ 
ronmental  effects  of  dredged  material  management  alternatives  -  a  technical  framework,”  EPA-842-B-92-008, 
Washington,  DC. 

Whelan,  G.,  J.  P.  McDonald,  J.  P.,  and  Sato,  C.  (1996).  “Multimedia  Environmental  Pollutant  Assessment  System 
(MEPAS):  Groundwater  pathway  formulations,”  PNNL-10907/UC-630,  Pacific  Northwest  National  Labora¬ 
tory,  Richland,  WA. 

Yeh,  G.  T.  (1981).  “AT123D:  Analytical  transient  one-,  two-,  and  three-dimensional  simulation  of  waters  transport 
in  the  aquifer  system,”  ORNL-5602,  Environmental  Sciences  Division,  Oak  Ridge  National  Laboratory,  Oak 
Ridge,  TN. 


NOTE:  The  contents  of  this  technical  note  are  not  to  be  used  for  advertising,  publication,  or 
promotional  purposes.  Citation  of  trade  names  does  not  constitute  an  official  endorsement  or 

approval  of  the  use  of  such  products. 


